Mode identification of high-quality-factor single-defect nanocavities 
in quantum dot-embedded photonic crystals 
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We investigate the quality (Q) factor and tlie mode dispersion of single-defect nanocavities based 
on a triangular-lattice GaAs photonic-crystal (PC) membrane, which contain InAs quantum dots 
(QDs) as a broadband emitter. To obtain a high Q factor for the dipole mode, we modulate the 
radii and positions of the air holes surrounding the nanocavity while keeping six-fold symmetry. 
A maximum Q of 17 000 is experimentally demonstrated with a mode volume of = 0.39(A/n)^. 
We obtain a Q/V of 44 000(n/A)'^, one of the highest values ever reported with QD-embedded PC 
nanocavities. We also observe ten cavity modes within the first photonic bandgap for the modulated 
structure. Their dispersion and polarization properties agree well with the numerical results. 



I. INTRODUCTION 

Photonic crystal (PC) nanocavities can confine light 
within a very small mode volume (V^) of less than a 
cubic wavelength with a high quality (Q) factor.^ Be- 
cause the light-matter interaction is enhanced in such 
structures, much attention has been paid to investigat- 
ing solid-state cavity quantum electrodynamics in PC 
nanocavities containing a semiconductor quantum dot 
(QD) as a quantum emitter intended for quantum infor- 
mation processing. Depending on the coupling strength 
between a cavity mode and a quantum emitter, physics 
of the system can be categorized into either a strong or 
weak coupling regime. In the strong coupling regime, the 
cavity and emitter coherently exchange energy back and 
forth (Rabi oscillation) j^i^i^ and this system can be uti- 
lized as a qubit. In the weak coupling regime, the spon- 
taneous emission rate of the emitter can be enhanced or 
suppressed compared to the rate obtained in a vacuum 
(Purcell effect) i^i^i This can be applied for not only effi- 
cient single photon source a^'^d^d^ but also polarization- 
entangled photon sources based on biexciton-exciton cas- 
cade emissions i^^d^d^ In the latter case, the nanocavity 
enhances the exciton decay rate to improve the fidelity of 
the entangled stateji^ii^ and its cavity mode should be 
polarization independent. 

Now we consider an optimum cavity structure for this 
purpose. There are two main types of PC lattice struc- 
tures: triangular and square. The first photonic bandgap 
(PEG) size of the triangular lattice is much larger than 
that of the square lattice. A wider PEG is preferable for 
a flexible design of high-Q nanocavities. Moreover, sev- 



eral cavity modes can exist within a wide PEG. In this 
case, the double resonance is possible with the proper 
cavity design, where one cavity mode is resonant with 
a ground exciton and another cavity mode with a first- 
excited exciton simultaneously. This enables an effec- 
tive resonant excitation of the exciton in a QD.^'' Several 
types of high-Q PC nanocavities based on the triangular 
lattice have been experimentally demonstrated, but six- 
fold symmetry is intentionally broken to obtain a high 
Q factor in most of these cavities J^ii^i^Sii^i These cav- 
ity modes are not polarization independent. The six- 
fold symmetry should be kept to apply cavities for the 
polarization-entangled photon sources^^^^^ 

In this paper, we present our study of the Q factor 
and the mode dispersion of single-defect PC membrane 
nanocavities made in a triangular lattice of air holes, 
which contain QDs as a broadband emitter. To obtain 
the high Q factor, we modulate the air holes surround- 
ing the cavity while keeping the six-fold symmetry. In 
sec. ini we present a numerical analysis of the Q factor 
and the cavity mode dispersion that we calculated for 
the modulated nanocavities. In sec. Illli the fabrication 
process of QD-embedded PC nanocavities and the exper- 
imental setup are explained. In sec. lIV Al we present our 
measurement of the Q factor of the dipole mode and com- 
parison with the calculated one. In sec. lIVEi we present 
our investigation of the characteristics of all the cavity 
modes within the first PEG. Ten modes are observed at 
maximum and their dispersion and polarization proper- 
ties are compared with the numerical results. Finally, a 
brief conclusion is provided. 
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FIG. 1: (a)-(d) Modified HI nanocavity structures, (e)-(h) Calculated Q factor as a function of modulation parameters of 
Ar/a and/or s/a for the structures of (a)-(d), respectively. The 3D-FDTD method was used with geometrical parameters of 
r/a = 0.31 and d/a — 0.71. A maximum Q of 33 000 was obtained at Ar/a = s/a = 0.09 as shown in (g). 



II. NUMERICAL ANALYSES 
A. Quality factor of dipole modes 

The nanocavity investigated in this work is a so- 
called HI nanocavity based on a triangular lattice of air 
holes formed in a thin GaAs membrane with one hole 
removed. ^•^ To obtain a high Q factor for the dipole mode, 
we modulated the radii and positions of the six nearest- 
neighbor holes to the cavity while keeping the six-fold 
symmetry."2iiSi2£i2J Though several ways can be used to 
perform the modulation, we simply changed the radii of 
the six air holes and/or their positions. Figures [Ha)-(d) 
show schematic diagrams of four types of modulated HI 
nanocavities. In these figures, a is a lattice constant and r 
(r') is a regular (modulated) hole radius. In Fig.[T](a), the 
six holes are shifted from the original positions (shown 
with dashed lines) to outside the cavity by s along the 
lines of symmetry with those radii fixed to r' = r. In Fig. 
[TJb), the radii of the six holes are reduced by Ar — r — r' 
without any position shifts. In Figs. [TJc) and (d), both 
the radii and position of the six holes are modified under 
the conditions of s = Ar for (c) and s = 2Ar for (d). 

The three-dimensional (3D) finite-diiference time- 
domain (FDTD) method was used to calculate the Q 
factors <^ The parameters included r = 0.31a, a slab 
thickness of d = 0.71a, and a refractive index of n = 3.4 
for GaAs. The size of the calculation domain was set 
to 31a X 16-\/3a in the slab plane and to 6a in the ver- 
tical direction to the slab, and the cavity was put at 
its center. The grid size was set to a/16. The calcula- 
tion domain was surrounded by the perfect matched layer 
(PML) based on Mur's second-order absorbing interface 
condition. The total energy stored in the cavity U(t) is 
expressed as follows: 

Uit)^Uocxp(-^t), (1) 



where ujq is the angular frequency of the cavity mode 
derived from the Fourier transform of the electric field, 
E{t). As can be seen in Eq. ([T]), the Q factor can be 
estimated from the slope of the exponential decay of the 
total energy within the cavity. Figures [Ue)-(h) show 
calculated Q factors as functions of modulation param- 
eters corresponding to Figs, [lla)-(d). As r' is reduced 
and/or s is increased, the Q factor increases drastically 
and then decreases. A maximum Q of 33 000 was ob- 
tained at Ar/a = s/a — 0.09 as shown in Fig. [Ijg). The 
mode volume was estimated to be F = 0.39(A/n)3 for 
the maximum-Q structure. 



B. Cavity mode dispersion and field distribution 

As mentioned earlier, the six holes were reduced and 
shifted to outside the cavity to obtain a high Q fac- 
tor for the dipole mode. Such modulations slightly en- 
large the cavity size, and accordingly, the energy of the 
cavity modes decreases. Figure shows the calculated 
cavity-mode frequencies as a function of Ar/a (= s/a) 
[see Fig. (He)] within the first PBG. The 3D-FDTD 
method was used again including geometrical parameters 
of r/a = 0.31 and d/a = 0.71. The solid (dashed) lines 
correspond to doubly degenerate (non-degenerate) cav- 
ity modes. Without modulation (Ar/a = 0), only the 
degenerate dipole modes exist within the PBG. As the 
radii decrease, several modes that originally exist in the 
second band without modulation appear within the PBG 
and move toward the first band. Ten modes including 
the degeneracy exist within the PBG at Ar/a > 0.075. 
A similar mode dispersion was reported by Park et al. 
theoretically.^ With the FDTD method, the energy of 
the cavity modes can be estimated from the Fourier 
transform of E{t). The degeneracy of the modes can 
be checked by introducing parity conditions. The plane- 
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FIG. 2: Normalized cavity-mode frequencies as a function of 
Ar/a (= s/a) calculated using the 3D-FDTD method. The 
solid (dashed) lines correspond to doubly degenerate (non- 
degenerate) cavity modes. The ten modes are referred to as 
modes A- J in order of frequency. 
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FIG. 3: Calculated electric field distributions [|i5|/max(|i5|)] 
in the slab center for ten cavity modes with Ar/a = 0.11. 
The modes in (a)-(j) correspond to modes A~J in Fig. [21 
respectively. 



wave expansion (PWE) method was also used to recon- 
firm it^i^ 

Figure [3] shows the calculated electric field distri- 
butions [|i?|/max(|i?|)] in the slab center for ten cav- 
ity modes with Ar/a = 0.11. Based on their mode 
shapes, they are caUed degenerate dipole [Fig. [31^a) and 
(b)], hexapole [Fig. [3jc)], degenerate quadrupole [Fig. 
[3]^d) and (e)], monopole [Fig. [3]^f)], degenerate second- 
quadrupole [Fig. OJg) and (h)] , and degenerate second- 
dipole modes [Fig. [3{i) and (j)], respectively, from the 
bottom to the top in the band diagram. For simplic- 
ity, the modes in Fig. [3Ua)-(j) are called modes A-J. 
In the real case, however, the degeneracy breaks due to 
small local fabrication error, which will be discussed in 



FIG. 4: Calculated electric field directions of (a) mode A, (b) 
B, (c) C, and (d) F. Modes A and B can be utilized for 
manipulating polarization-entangled states. Although non- 
degenerate modes of C and F have no preferred direction in 
the far fields, they cannot be applied for this purpose. 



detail later. To create a strong interaction between a 
cavity mode and a single QD, we should locate the QD 
at the antinode of the field. Since only a few groups 
can use a site-controlled QD,^ most use self-assembled 
or monolayer-fluctuated QDs whose positions cannot be 
controlled. Once a site-control technique is available, 
cavity modes that have a low number of nodes are desir- 
able because they can withstand a position error better 
than those with a lot of nodes. From this point of view, 
modes A, B, and F are preferable. 

Figure Hlja) and (b) show the electric field directions 
of modes A and B, respectively. These modes are ideally 
degenerate, and both modes have a strong electric field 
at the center of the cavity. Not only the far field directly 
above the cavity but also the local field at the center has 
no preferred direction. A single QD should be located at 
this position to manipulate polarization-entangled states. 
On the other hand, mode F is non-degenerate. This 
mode radially oscillates in the slab-plane with respect to 
the cavity center. Although the far field directly above 
the cavity is unpolarized, the local field at the antinode 
has a preferred direction, as shown in Fig. Iljd). There- 
fore, mode F can be applied for single photon sources 
and a strong coupling system but not for polarization- 
entangled photon sources. Mode C, whose electric field 
direction is shown in Fig.HJc), is also non-degenerate and 
the situation is the same as with mode F. 

Four modes are at the slightly higher energy side of 
mode F. For example, for a = 300 nm, the energy (wave- 
length) of mode F is 1.336 eV (924 nm), and the energies 
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FIG. 5: Specifications of fabricated QD-embedded PC 
nanocavity. (a) AFM image of the surface of a QD layer. 

(b) PL spectrum of the ensemble QDs at room temperature. 

(c) SEM image of the PC structure. The corresponding re- 
ciprocal lattice space is also shown. 

of mode G and / are 34 and 53 meV larger than that of 
mode F. These values are similar to typical energy sep- 
arations between the ground and first-excited excitons 
in self-assembled InAs/GaAs quantum dots.— Modes F 
and G (or /) may be simultaneously resonant with the 
ground and first-excited excitons, respectively. When a 
pump light is tuned to the first-excited exciton, it is en- 
hanced by the latter cavity mode and an electron-hole 
pair is effectively created only in the specific QD.^^ This 
resonant excitation suppresses the generation of excess 
carriers that degrade the quality of single and entangled 
photon sources. The ground exciton also interacts with 
the former cavity mode, which results in a Purcell effect 
or strong coupling. 



lective wet-etching techniquesi^ A reference QD sample 
grown with the same condition but without layers above 
and air holes was also prepared. Figure [^a) shows an 
atomic force microscope (AFM) image of the surface of 
the reference sample. The density of the QDs was esti- 
mated to be about 3 x 10^° cm^^ from the AFM image. 
The average height and diameter of the QDs were also 
estimated to be 4 and 30 nm, respectively. Figure (Hb) 
shows the photoluminescence (PL) spectrum of the fab- 
ricated QD ensemble at room temperature. We see not 
only the first peak (A ~f300 nm), which corresponds to 
the ground exciton transitions, but also the second peak 
(A ~I200 nm), which corresponds to the first-excited ex- 
citon transitions. Figure EI^c) shows a scanning electron 
microscope (SEM) image of a typical HI nanocavity. A 
number of nanocavities were prepared by changing struc- 
tural parameters, a regular air hole radius, r, and a mod- 
ulated air hole radius, r', while the lattice constant was 
fixed to a = 350 nm. The HI nanocavity was surrounded 
by 15 periods of air holes for good in-plane optical con- 
finement. 

The samples were set in a conduction-type microscope 
cryostat that cooled them down to 4.3 K. The cryostat 
was mounted on a piezo-actuated translation stage for 
precise sample positioningi^ The samples were optically 
excited by a He-Ne laser (A — 633 nm). The pump beam 
was focused by a microscope objective lens (numerical 
aperture of 0.42) to a spot size of 4 /im on the sam- 
ple. The pump power was set to ~3 kW/cm^. The 
PL was collected by the same objective lens and spec- 
trally resolved by a single-grating monochromator (focal 
length of 640 mm) equipped with a liquid nitrogen-cooled 
InGaAs multichannel detector (MCD, Jobin-Yvon IGA- 
512). The wavelength resolution of the system was 35 
pm at A = I3I0 nm with a I200-groove/mm grating. 
This was limited by the 50- /zm pixel width of the MCD. 
To obtain a higher wavelength resolution, we used a liq- 
uid nitrogen-cooled photomultiplier (PMT, Hamamatsu 
Photonics R5509-42) instead of the MCD, and both the 
entrance and exit slit widths of the monochromator were 
set to 20 /im. In the latter case, the wavelength resolu- 
tion of the system was 19 pm. 



III. EXPERIMENTS 

Based on the calculation results, we fabricated HI 
nanocavities with the modulation shown in Fig. [Ijc). 
They contain the QD ensemble which functions as a 
broadband emitter to probe the cavity modes. The sam- 
ples were grown by molecular-beam epitaxy on a (001) 
GaAs substrate. A 248-nm-thick GaAs membrane was 
grown on a 2-/im-thick Al0.gGa0.4As sacrificial layer on 
the substrate. A single layer of InAs self-assembled QDs 
were embedded in the middle of the GaAs membrane 
layer. The strain-relaxation layer was made on the QD 
layer to obtain longer wavelength emissions. After the 
growth of the membrane, HI nanocavities were made us- 
ing electron-beam (EB) lithography, dry etching, and se- 



IV. RESULTS AND DISCUSSION 

A. Quality factor measurement of dipole modes 

The Q factor of the dipole modes was extracted from 
the Lorenzian fit to the PL peaks measured at room tem- 
perature. Because the PL spectrum is the Fourier trans- 
form of Eq. ([1]), Q is equal to A/ A A, where AA is the full 
width at half maximum (FWHM) of the PL spectrum. 
Figure [H] shows the obtained Q factor as a function of 
the modulation parameter. The regular air-hole radius 
was r = 114 nm (r/a = 0.33). The modulation parame- 
ters of r' and s were estimated from the SEM image of 
the nanocavities. The PL-peak wavelengths were ^--^1300 
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FIG. 6: Measured Q factor as a function of the modulation 
parameter of Ar/a. A maximum Q of 17 000 was obtained at 
Ar/a — 0.11. The inset shows a spectrum of the maximum-Q 
nanocavity measured with the PMT. 
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FIG. 7: Typical spectrum of the dipole modes with Q > 
10 000. The small local fabrication error breaks the degen- 
eracy of modes A and B to make a 0.7-nm split. Inset: po- 
lar plot of the polarization dependence of short (open squares 
with line) and long (filled circles with line) wavelength modes. 
The geometry is the same as that shown in Fig. \^c) . 



nm, which correspond to the PL-center wavelengths of 
the QD ensemble. Five nominaUy identical cavities for 
each modulation parameter were investigated. As Ar/a 
increases, the Q factor rapidly increases and it becomes 
more than 13 000 at Ar/a = 0.10-0.11. The change in 
Q factor agrees well with the calculated one shown in 
Fig. [TJg). The cavity wavelength redshifts as Ar/a in- 
creases because the effective cavity size increases. This 
change is plotted in Fig. [HI The maximum Q of 17 000 
was obtained for one of the cavities with Ar/a = 0.11. 
The PL spectrum for this cavity taken with the PMT 
is shown in the inset of Fig. [6l The calculated mode 
volume, which is much less sensitive to the modulation 
parameter compared to the Q factor, is = 0.39(A/n)'^ 
for this cavity. The obtained Q/V, which is the figure 
of merit in the weak coupling regime, is 44 000(n/A)^. 
This is one of the highest values ever reported with QD- 
embedded PC nanocavities i^'^i^^'^^ The expected Purcell 
factor is Fp = 3Q{\/n)^ /Att^V = 3300 when a QD is 
located at the cavity center and when the cavity mode 
is right resonant with the QD exciton. The measured 
maximum-Q is lower than the calculated Q of 33 000 by 
factor of 2, and we think that the measured Q is limited 
by the QD absorption and the fabrication error through 
the EB lithography and dry etching process. 

Figure[7]shows a typical spectrum of dipole modes with 
Q > 10 000. The dipole modes A and B should be ide- 
ally degenerate. However, the degeneracy breaks due to 
the local imperfection in the fabricated struct urei^i^i 
The inset in Fig. [7] shows the polarization dependence 
of short (open squares with line) and long (filled circles 
with line) wavelength modes. Its geometry is the same 
as that shown in Fig. [5fc). They were linearly polarized: 
one is parallel to the F-K direction and the other to the 
F-M direction. The typical split was 6 = 1-2 nm. To 
analyze the effect of the small local fabrication error, we 



carried out an additional modulation to one of the six 
nearest-neighbor holes in the FDTD simulation. When 
the radius of the hole was reduced by 1 nm, the cavity 
modes split by 0.7 nm {5/X = 0.05%). Controlling the 
hole size with less than 1-nm error is difficult even with 
a state-of-the-art fabrication technique. Therefore, the 
degeneracy can be kept for low-Q (< 1000) cavities, but 
the mode split is unavoidable for high-Q (>10 000) cav- 
ities. A polarization-dependent cavity-frequency trim- 
ming technique, such as surface nano-oxidization by the 
AFMf2^ is required to compensate for it. 



B. Cavity mode identification within the photonic 
bandgap 

Figure [H shows that ten cavity modes theoretically ex- 
ist within the first PBG at maximum. The mode iden- 
tification of the fabricated structures was done by the 
PL measurement at 4 K. The wavelengths of the cavity 
modes and QD ensemble emissions were shifted to shorter 
wavelengths by about 20 and 90 nm compared to those 
at room temperature, respectively. 

Figure EJa) shows a wide-range PL spectrum of the 
high-Q sample with r/a ~ 0.31 and Ar/a = 0.11. The 
wavelength resolution of this spectrum is 0.1 nm. Ten 
sharp peaks were observed in broadband emissions of the 
QD ensemble. Though modes A and B are far detuned 
from the PL center of the QD ensemble, other emitters 
such as impurities are likely to contribute to the emis- 
sions. The calculated spectra obtained using the FDTD 
and PWE methods are shown in Fig. [HJb) and (c), re- 
spectively. In these figures, bold (thin) lines represent 
doubly degenerate (non-degenerate) modes. The mea- 
sured wavelengths agree well with the calculated ones. 
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FIG. 8: (a) Wide-range PL spectrum of high-Q nanocavity with Ar/a = 0.11 measured at 4 K. (b) and (c) Calculated spectra 
obtained using the FDTD and PWE methods respectively, where bold (thin) lines represent degenerate (non-degenerate) modes. 
The horizontal axis in (c) is shifted by —25 nm compared to the others in (a) and (b). (d)-(i) Polar plots of polarization 
dependence of modes A-J in (a). 



The horizontal axis in Fig. M^c) is shifted by —25 nm to 
enable the results to be compared easilyi^^ This is the first 
demonstration wherein ten cavity modes were clearly ob- 
tained for a single PC nanocavity. Conventionally, mode 
dispersion has been measured by plotting normalized fre- 
quencies taken from several samples having different ge- 
ometrical parameters. Dalacu et al. observed six cavity 
modes for an HI nanocavity, however, they did not ob- 
serve modes G-J.^^ The measured split for the ideally 
degenerate modes were 1.7 nm (modes A and B), 5.1 
nm [D and E'), 4.2 nm [G and iJ), and 3.9 nm (/ and 
J) . Because the positions of the antinodes in the higher 
modes are close to the modulated air holes [See Fig. [3], 
the local fabrication error affects the higher modes more 
than the fundamental modes. This resulted in a larger 
split of the higher modes compared to the fundamental 
modes of A and B. 

Figures [Sl[d)-(i) show the polarization dependence of 
the ten modes. In Fig.[5]Jg), the monopole mode F does 
not have a clear preferred polarization direction as was 



theoretically predicted. The non-circular shape of the 
data is mainly attributed to the polarization-dependent 
optical components. The hexapole mode C is linearly 
polarized along the F-K direction as shown in Fig. [SJe), 
though this mode should be unpolarized theoretically. 
This can be explained as follows: As can be seen in Fig. 
Sfc), this mode consists of three identical components 
that oscillate along the F-K direction. The three com- 
ponents combined equally, and the mode does not have a 
preferred direction when six-fold symmetry is held. Once 
the symmetry breaks, one of the three components is 
dominant, and the mode is polarized parallel to one of 
the F-K directions. Modes A and B, G and //, and / 
and J are pairs of orthogonal polarization [Figs.IDJd), (h), 
and (i)], respectively, whereas modes D and E are not or- 
thogonal but their polarization directions differ by ~45° 
[Fig. IHI^f)].'^^ This non-orthogonal polarization property 
was observed in our similar structures and elsewhere 

Figure [5] shows the measured cavity mode wavelengths 
as a function of the modulation parameter. The calcu- 
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lated wavelengths are also shovifn with solid and dashed 
Hnes. The measured ones agree well with the theoretical 
ones. Though the simulation predicted that ten cavity 
modes exist at Ar/a > 0.075, we experimentally ob- 
served ten cavity modes only at Ar/a > 0.10. This is 
mainly because the bandedge frequency of the second 
band of the fabricated structure is lower than the calcu- 
lated one. This also affects the Q factor of modes G- 



J. Because the cavity mode is closer to the bandedge, 
the in-plane optical confinement decreases, as does the 
Q factor. The measured Q of modes G-J was ^^1500 at 
Ar/a = 0.11, which is more than 10 times lower than the 
calculated one. 



V. CONCLUSIONS 

We have theoretically and experimentally investigated 
the Q factor and mode dispersion of QD-embedded HI 
nanocavities. By modulating the six air holes surround- 
ing the nanocavity while keeping six-fold symmetry, we 
experimentally obtained a maximum Q of 17 000 for the 
dipole mode with a very small V = 0.39(A/n)'^. These 
values correspond to a Purcell factor of Fp = 3300 under 
an optimum condition. The measured Q was compatible 
with a calculated Q of 33 000. Though the dipole mode 
was designed to be doubly degenerate, the small fabri- 
cation imperfection slightly broke the symmetry, which 
caused a mode split of about 0.1%. Ten cavity modes 
were observed in the single HI nanocavity within the 
first PBG, and their dispersion and polarization prop- 
erties agreed well with those simulated. The measured 
energy difference between some higher modes was almost 
the same as that between the ground and first-excited ex- 
citons in a QD, which would meet the doubly resonant 
condition. Once the cavity modes are precisely tuned to 
the exciton levels in a QD, efhcient non-classical photon 
sources can be achieved for quantum information pro- 
cessing. 
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